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Abstract 

T his report discusses different approaches to stable H°° controller design applied to the 
problem of augmenting the longitudinal dynamics of an aircraft. Stability of the H controller 
is investigated by analyzing the effects of changes in the performance index weights, and 
modifications in the measured outputs. The existence of a stable suboptimal controller is 
also investigated. It is shown that this is equivalent to finding a stable controller, whose 
infinity norm is less than a specified bound, for an unstable plant which is determined from 
parametrization of all H°° controllers. Examples are given for a gust alleviation and a command 
tracking problems. 


‘This work was performed while the first author was a NASA/OAI Summer Faculty Fellow at NASA Lewis 
Research Center. 




1 Introduction 


The problem of finding a stable feedback controller is known as strong stabilization, and it has 
been widely studied in the literature, see e.g. [1]. One of the motivations behind this problem is 
that, when the plant is stable, unstable controllers are not tolerant to faults in measurements. 
An example of this is when a feedback path is broken, such a controller may lead to an 
unbounded response for a bounded reference input. Also, from a real-time implementation 
point of view, it may be undesirable to have an unstable feedback controller, [2], [3]. 

A necessary and sufficient condition for the existence of a strongly stabilizing controller 
is the parity interlacing property (p.i.p.) [4]. A plant satisfies p.i.p. if the number of poles 
between any pairs of distinct blocking zeros on the positive real axis (including + 00 ) is even. 
There are procedures for constructing stable controllers which stabilize a given plant, [5, 1, 4]. 
However, there is no simple parametrization for the set of all strongly stabilizing controllers. 
Usually, closed loop stability is the first design requirement. But one is also interested in 
achieving some kind of robustness and performance level in the controller design. This require- 
ment can be achieved (with a certain degree of conservatism) by using 7i°° control techniques. 
Therefore, finding strongly stabilizing controllers is an important research problem, which 
is the subject of this report. It should also be mentioned that some promising results appear 
in [6] (see also references therein) on the version of this problem. The effects of weight 
selection on the stability of the optimal W°° controller for SISO plants have also been studied 



The plant considered in this study is a linear model of the longitudinal dynamics of an 
experimental F-15 aircraft with pitch vectoring nozzles. For this plant, two 7f°° control de- 
sign problems are defined. The first one deals with gust alleviation, and the second one is 
a command tracking problem. The results of [8], [9] and [10] are used in order to obtain a 
parametrization of all suboptimal controllers. The underlying operator theoretical results 
for this parametrization can be found in [11, 12], see also [13] for more details. Most com- 
mercially available softwares (e.g. robust control modules of MAT LAB and MATRIXx ) 
generate the so-called “central controller” of [9], (see Section 2.2 for the precise definition of 
the central controller) . Here, the effect of structural changes in the plant (e.g. adding one 
more output for feedback) on the stability of the central controller is studied first. Then, the 
effects of scaling the performance index weighting functions, and the effect of increasing (or 
decreasing) the 7 i°° suboptimal peformance level 7, on the stability of the central controller 
are studied. Finally, the parametrization of all suboptimal H°° controllers is studied, and 
several different methods of finding a stable controller in this parametrization are discussed. 
This is a significant problem, because for a given admissible 7 the central controller may be 
unstable, but there may be a stable controller in the set of all controllers which achieves the 
same performance level. 

The rest of this paper is organized as follows. In the next section a background on the 
standard H°° control problem is given, along with the formulae for the central controller and 
the parametrization of all suboptimal controllers. The plant model considered in this paper is 
also described in the next section. Section 3 is devoted to the studies on the stability of the 
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central controller. In Section 4, stable suboptimal H°° controllers are investigated for the gust 
alleviation and tracking examples. For both of these examples it is shown that there exists 
a stable suboptimal controller, while the central controller is unstable. Finally, concluding 
remarks are made in Section 5. 

2 Background on 7 i°° control 

2.1 Standard problem set-up 

The so-called “standard H°° control problem” deals with the system shown in Figure 1. The 
system equations are assumed to be given by the following 


x(t) = 

Ax(t) + Biw(t) + B< 2 u(t) 

(i) 

z(t) = 

C\x(t) + D u w(t) + Di 2 u(t) 

(2) 

y(t) = 

+ ^21^(0 + ^22^(0 

(3) 


where x represents combined states of the system, and components of w are the exogenous 
signals (reference inputs, disturbances, measurement noises), components of u are the control 
inputs, components of y are the measured signals, and components of z are internal signals to 
be controlled. The optimal 7i°° problem is to find a feedback controller K (whose input is y and 
output u) so that the closed loop system is stable, and the worst energy amplification from w to 
z is minimized. This problem is equivalent to finding a stabilizing controller which minimizes 
Halloo, where T zw (s) is the closed loop transfer function from w to z. The suboptimal 
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7i°° control problem is to find a stabilizing controller so that HTzwIioo < 7 , for a specified 
performance level 7 . 

Usually, an 7f°° control problem is first transformed to the above standard form. Then, 
the controller is obtained from the matrices A, B \ , B 2 , C\, C 2 , Du, Du, D 21 , D 221 using 
algebraic Riccati equation solvers. The controller formulae are given in Section 2.2. 


2.1.1 Aircraft Model 


In order to demonstrate how one sets up an control problem, two examples will be con- 
sidered. Both of these examples involve a nominal plant, which is a linear model for the 
longitudinal dynamics of an aircraft. The plant is decribed by 


x p (t ) = A p x p (t) + BpU v (t) 

y?(0 = CpXp(t) 


where x p are the states, u p denotes the command input and y p denotes the plant output, and 
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c p = 


where C a = 57.2958[ 0 10 0 ], C qp = 57.2958[ 0010 ] 


L'-'qp J 


The first component of y p (t) is the angle of attack a(t) (in degrees) and its second component 
is the pitch angular rate q p (t) (in degrees/sec). The command input u p is the nozzle pitch 
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vector angle (in degrees). The components of the states x p of the airframe are forward velocity 
(in ft/sec), angle of attack (in radians), pitch rate (in rads/sec), and pitch attitude 9 (in 
rads), respectively. In this system the command signal u p is the output of an actuator whose 
dynamical behavior is decribed by 

x a (t) = —25x a (t) + 25u(t) 

Uji(0 = ®a(0 i 

where u(t ) denotes the command signal to be generated by the feedback controller. 

The poles of the transfer function T Upt y p , from u p to y p , are 

[+0.634, -0.8647, -0.2206 ± j0.1274] , 

and the zeros of the first component T UpiQ and second component T Up<qp , of T UptVp are 

[-39.6, -0.1382 + jO.1726] and [-0.436556, 0.000, +0.002912] 

respectively. Note that this plant satisfies the parity interlacing property, but the transfer 
function from u p to q p does not. In other words, there exists a strongly stabilizing controller 
for this plant. But, if the first output a(t ) is not used as a control feedback, then for the 
resulting plant all stabilizing controllers are unstable. 

2.1.2 Tracking Problem Definition 

Now the H°° control problem, associated with tracking of a-command signal, can be defined in 
terms of the system shown in Figure 2, where Ws(s) is the sensitivity weight, and VFx(s) is the 
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weight for the complementary sensitivity. Usually Ws is chosen as a low pass filter, and Wj 
is chosen as a high pass filter. This means that low frequency reference inputs are considered 
for tracking, and high frequency unmodelled dynamics are taken into account for robustness, 
see [5]. Accordingly following weights are chosen in this problem 

Ws(s) = C s (s - A a )~ l B s where A s = -1/335, B s = 1000/335, C 3 = 1, 

W T {s) = Ct($-A t )- 1 B t + D t where A t = - 1/.00278, B t = -99.999/.00278, C t = 1, 


and D t = 100, with R x = 1/3, R 2 = 1/3, R 3 = 0.2, R 4 = 0.02 Ki = 3, K 2 = 0.001. The 
scalars Ri,...,R 4 assign relative weights on the internal signals of interest and K 1 , K 2 scale 
the measurement noises. For example if the size of u p has to be made small then R 3 (weight 
on Up) is chosen relatively large. In Figure 2, represents the reference input (we want a to 
follow u>i) and K 2 w 2 represents a small amount of noise which is present in the measurement 

f V\jR\ 1 


of q p . If y := 


and u := 2bR 4 u, then the state equations, corresponding to the 


U 2 /a' 2 J 

generalized system from [iv u] to [z y], are in the form (1-3) with 
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Note that the controller generates u = K(s)y. But first the controller which generates u 

~K i Ol" 1 

K(s)v , can be found and then from K(s) one can determine K(s) = (25 J^)" 1 

L 0 il 2 J 

Numerically it is simpler to find Jf, because in this set-up and D 21 are normalized to have 


entries 1 and 0 only. See the formulae in [9] for a comparison. 


2.1.3 Gust Alleviation Problem 


Consider the plant described above with gust affecting the system dynamics as follows 

"wi (O' 


— ApgXpg(t} + Bg 
y V (t) = CpgX pg (t ) 


L^2(0 




where C pg = 


Apg — 


0ix3 C Q 
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, in other words y p = 
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n o oi 
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with V = 300ft/sec, L = 1670ft, and 0 < a < 10, (a moderate gust, a = 5, will be considered). 
The gust model in the above formulae is the Dry den model, which is described in [14]. The 
command signal u p is generated by an actuator whose input is u, which is to be generated by a 
feedback controller: u p = W^v., where W^s) = C a {$ — A a )~ l B a , with A a = —25, B a = 5 and 
C 0 = 5. The Ti 00 problem associated with this system is to minimize the effect of w\ , . . . , w 4 
on zi, . . .Z 4 as shown in Figure 3, where W 3 /K 1 and W 4 /K 2 represent measurement noises, z\ 
and Z 2 are scaled values of a and q p respectively, and 23 is the scaled control, and z 4 is a blend 
of the control and the control rate. 


For this problem define y := 


A'lJ/il 


, then the system equations which represent the 


LA' 2 y 2 J 

transfer function from [u>, u] to [ 2 , y] are in the form (1-3), with 
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C2 = 


0 


K\ C ag 

K^Cqpg . 


D 2 l = [02x2 12x2 ] 


and £>n = 0 4X 8, -D22 = O2 x 1 * Similar to the previous case, once a controller K , which generates 
u from y, is determined, the equivalent controller K , which generates the same control u from 


y , can be found by setting K(s) = K(s) 


K\ 0 ' 
L 0 K 2 . 


2.2 Controller Formulae 


The formulae of [8, 9, 10] for 7f°° controllers, which satisfy a certain specified performance level 
7, is given below. The problem formulation, in both tracking and gust alleviation, satisfies the 
structural assumptions of [9]. The controller expression is simpler if the formulation satisfies 
the structural assumptions of [8]: 


Assumptions 

Al. (A, B 2 ,C 2 ) is stabilizable and detectable, 
A2. Du — 0, D 22 = 0, 



' Bi - 


'O' 

A3. 


Dl 1 = 



.£21. 
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A4. 


[Cl Du] = [0 /]. 


The gust problem defined above satisfies these assumptions. For the tracking problem the 
orthogonality assumptions A3 and A4 are not satisfied, therefore for this example a slightly 
more complicated formula, [9], will be used. 

Now consider the problem formulation given in (1-3), and assume that they satisfy Al, 
. . ., A4, e.g. the gust alleviation problem is in this form. Then, set up the following Algebraic 
Riccati Equations (AREs) 

A t X +XA + Xi^rBtBf - B 2 Bf)X + C{Ci = 0 (4) 

7 

AY + YA t + Y(^Cf Ci - Cl C 2 )Y + BiBl = 0 . (5) 

7 

It has been shown that, [8, 9], there exists a stabilizing controller which makes HT^Uoo < 7 if 
and only if there exist unique solutions X and Y , for (4) and (5), such that 

(i) X and Y are positive semi-definite, 

(ii) eigenvalues of Ax := A + {^B\Bl — B 2 Bl)X and Ay := A + Y(-^ClCi — ClC 2 ) are 
in the open left half plane, and 

(iii) the largest eigenvalue of Ay is strictly less than j 2 . 

Moreover, if (i)-(iii) hold then all controllers, K, which stabilize the closed loop system and 
achieve HT^Hoo < 7 are parametrized as follows 

x c (t) = A c x c (t) + B ic y(t) + B 2c q{t) (6) 
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u(f) — Cic^-cCO -h ?(0 


( 7 ) 


r(t) = C 2c x c (t) + y(t) (8) 

where 


A c = 

A X - ZYC?C 2 , 

( 9 ) 

B lc = 

zrcf, 

(10) 

B 2c = 

zb 2 , 

(11) 

C u = 

-Bjx, 

(12) 

C 2c = 

-c 2 

(13) 


and q(t) is generated by a transfer matrix Q(s), whose input is r(f), with Q € V. 00 (that is 
Q must be stable) and HQHoo < 7. Obviously, there are infinitely many choices for Q(s) and 
hence there are infinitely many suboptimal controllers. Implementation of this controller is 
shown in Figure 4. In particular, one can choose Q(s) = 0, this gives the “central” controller. 
The issue of selecting Q will be discussed in Section 4. In [9] and [10] an alternative expression 
for suboptimal controllers is given. This expression is in the same form as (6-8), but its state 
space realization is slightly different 


Ac 

= Ay — B 2 B 2 X Z, 

(14) 

B\ c 

II 

a. 

(15) 

B 2c 

= B 2 , 

(16) 

Clc 

= -BjXZ, 

(17) 
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C*2c — ~C 2 Z. 


(18) 


The above realization (14-18) can be obtained via a state transformation on (9-13). 

It should be mentioned that for the tracking problem, the parametrization of all suboptimal 
controllers is exactly the same as (6-8), but the AREs associated with this problem, and hence 
computations of A e , B \ c , £20 C\ c and C20 are slightly more complicated, see [9] for all the 
details. 

Note that the central controller is stable if and only if A c is a stable matrix (i.e. all the 
eigenvalues of A c have strictly negative real parts). Whereas, stability of a controller which is 
obtained from a non-zero Q(s) depends on A e , l? 2c , C 2c and Q(s), to be discussed further in 
Section 4. 

3 Stability of the central controller 

In this section the effects of structural changes and weight scalings, on the stability of the 
central controller, will be discussed. 

3.1 Changes in the output variables 

Consider the tracking problem defined in Section 2.1.2. The central controller corresponding 
to 7 = 2 is unstable (for this system the optimal value of 7 is between 1.5 and 2). It will be 
shown that by changing one of the output variables it is possible to transform the problem to 


12 



a form where it is easy to analyze stability of the central controller. For this purpose replace 
the second output 2/2 = 9p + K 2 w 2 with y 2 — u p — a where u p = C a x a + K 2 w 2 . Here K 2 id 2 
represents a small amount of actuator noise, whereas in the previous case it was a measurement 
noise. Also for simplicity redefine 24 from R 4 U to 24 = u. Then, the generalized plant is as 
shown in Figure 5. 

As before, if ft := yi/K for i = 1,2, then the system equations are as in (1-3) where 



' A a 
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0 D t C a 

0 Ct 

c 2 = 

■ 0 

—C a /K\ 0 0- 
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. 0 
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0 . 





and D\z =[0 0 0 l] r , D 22 = /, Du = 0, D 2 2 = 0. It is easy to check that if below assump- 
tions (dffl-dfF3) are satisfied then the problem formulation is in the form of the disturbance 
feedforward (dff) problem defined in [8], 

dffl. C a B a is invertible, D t ± 0, K 2 ^ 0, 
dfF2. ( A P ,B P , C Q ) is stabilizable and detectable, 
dfI3. A a , A s , A t and {A v + B p C a ) are stable. 
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In this case the A-matrix of the central controller is 


A c := A - B X C 2 - B 2 BlX. 


For the plant shown in Figure 5, the above dff assumptions are satisfied. Indeed one can check 
that for A p , B p and C a given above the eigenvalues of (A p + B p C a ) are at —0.1221 ± j'0.1757 
and —0.2597 ± j’1.7419. It is also easy to verify that in this case A c is in the form 
•A a -B a BjX u * N N" 

0 Ap + B p C a 0 0 


A c = 


0 0 A, 0 

0 BtC a 0 At. 


where X n is the 1, 1 entry of X , and N represents an entry whose value is not important. By 
partitioning this matrix one can easily see that the eigenvalues of A c are the eigenvalues of its 
diagonal blocks. Since A p + B p C a , A 3 , A t , A a are stable and B a BjX u is a positive scalar, 
one concludes that the central controller in this case is stable. 


The purpose of this example was to show that by modifying one of the outputs it was 
possible to obtain a central controller whose A-matrix A c has a certain special structure for 
which the stability analysis is easy. Of course for more complicated examples it may be difficult 
to see which way the output should be modified to get a stable controller. Moreover, in some 
cases it may be impossible to modify the output physically. 
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3.2 Scaling the weights and adjusting 7 


Now the effects of weight scalings on the stability of the central controller will be studied for 
the gust alleviation example shown in Figure 3. Note that in this example one can think of K 1, 
K 2 , Ri, R 2 , R 3 , Ra and a as scaling parameters. As illustrated by Figures 6 and 9 through 11, 
changes in K\, R\, R 2 and R 3 do not significantly affect the stability of the central controller. 
These figures are generated with the nominal values R\ = R\ = R 2 — R 3 = 1> — 0-01 and 

a = 5. It is seen from Figures 7 and 8 that with these nominal values, if K 2 is larger than 0.1, 
the central controller is unstable, for both 7 = 25 and 7 = 30. In the above analyses, large 
enough 7 values are chosen in order to guarantee the existence of an H°° controller for the 
parameter variations considered in each plot shown in Figures 6 through 11. 

Now fix K\ = K 2 = R\ = R 2 = R3 = 1 and <7 = 5. For these values of the scaling 
parameters 

18 < 7 opt < 19 , 

and the eigenvalues of the central controller for two different values of 7 are given below. 


15 



mam 

7 = 25 

-0.16±;0.06 

—0.17 ijO.053 

-0.18 ± jO.0057 

-0.18 ± j'0.00567 

+ 0.6 

+ 0.3 

-10.5 

-5.5 

-16.8 

-17.1 

-34.3 

-34.9 


This example shows that by increasing the value of 7 it is possible to move the unstable pole 
slightly towards the left half plane. But the central controller remained unstable even for 7 as 
large as 80. In fact, it is known that [8], as 7 — ► 00 the central Ti°° controller converges to the 
optimal 7 i? controller, which does not have to be stable. 

In the tracking problem, the general rule of thumb is that if one is willing to give up the 
performance to gain more robustness, then one should scale the sensitivity weight Ws by 1 /k 
and complementary sensitivity weight by A: > 1. Then, for k sufficiently large the central 
controller is expected to be stable, see [7] for a delay system example. To illustrate this point, 
now consider the generalized plant of Figure 12, which is a slightly different version of the 
tracking problem defined in Section 2.1.2, (there are two more outputs for feedback, a and 6, 
with slight additive noise). For this problem when k = 1 the optimal value of 7 is between 1.5 
and 1.75, and for 7 = 1.75 the central controller is unstable. By increasing Wj and decreasing 
Ws by a factor of k > 2 one gets a stable central controller for 7 = 2, as illustrated by 
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Figure 13. However, this decreases overall performance, and increases response to noise. 


4 Existence of a stable suboptimal controller 


The previous section was devoted to studies on the stability of the central controller for a given 
performance level 7. Sometimes the central controller may be unstable, but there may exist 
another controller which is stable, and achieves the same performance level 7. In Sections 4.1 
and 4.2 this situation will be illustrated with the gust alleviation and tracking examples. But 
first some possible ways to find such a controller are described. 

Recall that, for a given 7 > 7 opU all suboptimal controllers are parametrized by (6-8) where 
q = Q(s)r, with Q e H°° and ||Q||oo < 7- That is Q is the free parameter, but it has to be 
stable, and its norm has to be bounded by 7 . Suppose that Q has a state space realization 

Xg(t) = AgXg(t) + Bgr(t) 
q(t) = CgXg(t) + D q r(t) . 


Then, inserting q{t) into (6) and (7), and using r(t) given by (8), one gets the following 

realization for the controller (which is a transfer function from y to it) 

Ac + B2c-DgC2c B2cCg 


’i c (ty 


.Xg{t ) . 



B q C2c 


U (t) = [ C\c 4 " DqC2c Cq] 



■*c(*r 


B\ c + B2 c Dg 



+ 


J 



- Bg . 


y(t) 


X c {t) 

U 9 (0J 


+ D q y{t) . 


(19) 

( 20 ) 
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Note that the A-matrix of this controller is 


Ak := 


A c + B 2c D q C 2c 
B q C 2c 


B 2c G q 

A q . 


( 21 ) 


Hence, there exists a stable suboptimal H°° controller for a given performance level 7 if and only 
if there exists a “subcontroller” Q{s), stabilizing the “subplant” G c (s) := B 2c (sl - A c )~ l C 2c , 
such that Q € ft 00 and ||<?||oo < 7- Thus, testing whether there exists a stable suboptimal 
H°° controller amounts to checking if G c is strongly stabilizable by a controller whose norm is 
bounded by 7. In other words, it is necessary (but not sufficient) that G c satisfy the parity 
interlacing property. If G c satisfies p.i.p., then the next question is: can we find a norm 
bounded stabilizing controller? To the authors’ knowledge there have not been any published 
results on this problem. 


One possible way to approach this problem is to apply any stable controller design procedure 
(see e.g. [1, 5 , 4]) to G c {s), and check whether this controller’s oo-norm is bounded by 7. 
Obviously, there is no guarantee that this method will work. A rigorous approach would be 
to study different stable controller design algorithms, and try to determine if it is possible to 
modify them so that the oo-norm of the resulting controller is bounded by 7. 


Another approach, which is more direct, would be to study the structure of all controllers 


stabilizing G c , and determine ways to choose the free parameter so that the controller is stable 
and norm bounded. To illustrate this method assume that Gc is SISO, then all stabilizing 
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controllers are in the form [1] 


0 _ Xg + DgQq _ Dg(Xg + £gQq ) 

Yg — NgQq 1 — Ng{Xg + DqQq ) 

where G c — Ng/Dg and Ng, Dg, XgiYg € W 00 satisfy the Bezout equation 

N g X g + Y g Dg = 1 , 

and Qq € 7f°° is the free parameter. So, in order to have Q G H°° and ||Q||oo < 7> one should 
try to make ||X G + D g Qq\\oo as small as possible. However, by the Bezout equation X G is 
equal to 1 /Ng, at the zeros of Dg in the right half plane. So there is a limit how small this 
expression can be made. Because of time limitations the authors did not further investigate 
this method of trying to select an appropriate Qq € . It may or may not give satisfactory 

results. 

Note that, when the realization (19-20) is minimal, a dynamic Q(s) (i.e. Q(s) is non- 
constant) leads to a controller whose order is larger than the order of the generalized plant. 
This may be an undesirable situation. Therefore, one may want to restrict the search to non- 
dynamic Q's, in other words Q(s) — Dq which is a constant. In this case, Q is automatical!} 
in H°°, so the only restriction is on the norm, that is o max (D q ) < 7. Moreover, if the number 
of inputs y and outputs u, of G c , is small then a parameter search, in the space of the entries 
of D q , can be done to find a feasible D q . In the next section this approach will be illustrated 
on gust alleviation and tracking examples defined earlier. 
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4.1 Example: gust alleviation problem 


Consider the 7i°° control problem defined in Section 2.1.3, with the parameters Hi = R 2 = 
R 3 = K\ = Ji2 = 1 and a — 5. For this problem 18 < 7 opt < 19; and for 7 = 25 the central 
controller is unstable, with a right half plane pole near +0.3. Now suppose 7 = 25 fixed. Then, 
the problem is to determine if there exists D q = [D q 1 D q 2] with 

D] 1 + D 2 q2 < 25 2 (22) 

such that the eigenvalues of Ak = A c + BicDqCic are in the open left half plane. The region in 
the D 9 -plane which gives an affirmative answer to this question is the shaded area in Figure 14, 
which falls in the circle (22). Note that allowable D q 's are in the form 

D q = R [cos 6 sin 6} 

where 0 < R < 25 and -7r < 6 < tt. The plot of maxReA(AA'), as R varies between 20 and 
25, are as shown in Figure 15, for two cases of fixed 0 = —10° and 6 = —26°. Also, the plot of 
maxReA(AA'), as 6 varies between -40° and -28°, are given in Figure 16, for fixed R at 23.5 
and 24.9. 

Figures 15 and 16 show that 0 should be small and R should be large in order to make Ak 
stable. Besides the stability of Ak-, it is necessary to check the overall system performance in 
the frequency domain, as well as in the time domain, in order to determine the optimal value 
of D q . To illustrate the overall system performance allow to be set equal to D q — [24.75 0.00]. 
As expected the closed loop performance is within specified level 7, i.e. cr max {T zw {j u>)) is below 
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7 = 25, see Figure 17. For this system, the effects of a pulse gust (i.e. = w 2 (t) are equal 

to a unit pulse of duration 0.25sec) on z\ , . . . , z± are shown in Figure 18. 

4.2 Example: tracking problem 

For the tracking problem defined in Section 2.1.2, recall that 1.5 < 7 opt < 2, and for 7 — 2 
the central controller generated by the algorithm given in [9] is unstable. Now consider the 
problem of determining if there exists D q = [D q \ D q 2] such that Ak = A c + B 2c DqC 2 c has all 
its eigenvalues in the open left half plane, and 

D 2 ql + D 2 q2 < 2 2 . ( 23 ) 

The shaded region in Figure 19 which falls in the circle (23) indicates possible values of allowable 
D g which make A K stable. Note that compared to previous example, we don’t have much 
freedom in choosing D q , it has to be chosen close to D q = [-1.99 0.00] in order to get a stable 

Ak- 

In the remainder of this section the performances of the suboptimal 7i°° controller obatined 
by choosing <2(s) = D q = [-1.99 0.00], and the unstable central controller, i.e. the one with 
D q = [0 0], will be compared. Frequency domain closed loop performances are compared in 
Figure 20. As expected both controllers yield UT^Hoo < 2. 

It should also be mentioned that the observation noise A' 2 tu 2 is added to our original H°° 
control problem definition in order to satisfy rank conditions required in the solution procedure 
given in [9]. Also, the weights are added to the standard H°° control problem in order to obtain 
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a controller, which gives desired loop shapes; i.e. the weights are not physically present in the 
system. So, once a controller is obtained it can be implemented on the physical plant as shown 
in Figure 21. In this figure a c (t) represents the reference (command) input to the controller, 
which generates u. The signal u is an input to the actuator, whose output u p is the plant 
input. As before, plant outputs are a and q v . 

Time domain performances of the controllers obtained for 7 = 2, with D q = [-1.99 0.00] 
and with D q = [0 0], are shown in Figure 22. It is seen that the response a(t) is better 
for D q = [0 0] (i.e. better tracking of a negative unit step o c (0)» the command signals 
u p (t) are similar for both cases. Frequency domain performance plots are shown in Figures 23 
and 24, where individual Bode plots (from a c to u p (t) and a(t) respectively) are shown. As 
expected the main difference in the frequency domain characteristics are in the high frequency 
behavior. Note that when D q ^ 0 the controller is not strictly proper, so compared to D q — 0 
there is less high frequency roll-off. 

5 Conclusions 

In this research possible ways of obtaining stable controllers, achieving a certain pre-specified 
H°° performance level 7, are identified. First the central controller A c is made stable by 
scaling the weights and/or modifying the measured outputs. Secondly, existence of a stable 
suboptimal controller is investigated by adjusting the free parameter Q(s) which appears in 
the characterization of all suboptimal controllers. In the first method, since the structure of 
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the H°° problem is modified, the performance of the stable controller may be different than 
the original unstable controller. In the second method, it is expected that if 7 is close to 
its optimal value, then all suboptimal controllers are expected to be close to the central 
controller, see e.g. [15]. Therefore, if the central controller is unstable and if 7 is close to 
there may not exist a stable suboptimal controller. So, for the second approach to work 7 has 
to be “sufficiently larger” than 7 opt . 

The second approach leads to an interesting problem: given an unstable plant G c (s), find 
a strongly stabilizing controUer Q(s) such that ||<?||oo < 7- This problem is difficult to solve 
in complete generality. Therefore the authors have concentrated their efforts in demonstarting 
the existence of a constant norm bounded Q(s) — D q which stabilize (j c (s), for both gust 
alleviation and tracking problems. The approach here was simple brute-force. Since in this case 
the number of inputs and outputs, for the controller, was 2 and 1 respectively, the parameter 
space was two dimensional; and it was easy to check sufficiently large number of points to get 
an idea about the region in the D q - space that gives a solution. 

Recall that the controller A-matrix, when D q is used, is given by Ak = A c + B 2 cD q C 2 c, 
and from (9-13) A c = A x - ZYCjC 2 , B 2c = ZB 2 , and C 2c = -C 2 . Therefore, 

Ak — Ax — Z(Y C 2 + B 2 B q)C2i (^4) 

where Ax is stable. This structure might give an idea about how D q should be chosen for 
more complicated input/output definitions. For example, one may try to minimize the norm 
of ZiYCj + B 2 D q )C 2 - This is a relatively easy optimization problem, which can be solved 
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using existing commercially available software. However, there is no guarantee that after 
minimization all the eigenvalues of Ak will be in the open left half plane. 

The effects of weight scalings and modifications in the measured outputs, on stability of 
the central controller, are studied in the first part of this research. Recall that A c is in the form 
(9) or (14). More specifically A c = Ax - ZY Cj Ci or A c = Ay — BiB^XZ where Ay and 
Ax are stable matrices. If the measured outputs are modified, then Cj changes, and without 
changing Ax it changes the term ZYCJC 2 . Therefore, one might be able to design C 2 so that 
ZYCJC 2 is “small” compared to the eigenvalues of Ax, and by using Gershgorin’s theorem 
[16] it may be possible to guarantee a stable A c . Similarly, if control inputs can be redefined, 
i.e. j f? 2 can be modified, then one can try to make B 2 BJXZ “small” without changing Ay, 
and this may lead to a stable controller. It is obvious from (24) that simultaneously modifying 
C 2 and selecting D q which minimizes Z(Y Cj + B 2 D q )C 2 increases chances of getting a stable 
controller whose poles are close to the eigenvalues of Ax ■ 
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Figure 1: Standard Feedback System 
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Figure 2: Tracking Problem 
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Figure 3: Gust alleviation problem 
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Figure 5: Alternative tracking problem definition 
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Figure 7: The effect of Kj on the most unstable eigenvalue of A c ; for y = 25 
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Figure 8: The effect of Kj on the most unstable eigenvalue of A c ; for y = 30 
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Figure 1 1 : The effect of on the most unstabl 
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Figure 12: Tracking problem with more outputs, and scaled W 5 , Wj 
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Figure 13: max ReA(Ac) versus A;, for 7 = 2 
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Figure 16: Effect of 6 on the most unstable eigenvalue of Ak, for R = 23.5 and R = 24.9 
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Figure 17: Closed loop performance for the gust alleviation problem; 
“maxCfzwd 40 )) versus co, for D q = [24.75 0.00] 
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Figure 18: Time resposes to a unit pulse gust, for D q = [24.75 0.00] 
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Figure 19: Values of D q which make Ak stable 
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Figure 20: Closed loop performance for tracking problem; amaxCTzwO 4 * 1 )) 
versus <o, for Dq = [-1.99 0.00], and for Dq = [0 0] 
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Figure 21: Controller implementation on the physical plant 
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Figure 22: Time responses to a negative unit step r c (<); 

Dq = [—1.99 , 0] (solid line), Dq =[0,0] (dashed line) 
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Figure 23: Bode plot of the closed loop transfer function from a c to u p ; 
Dq = [-1 .99, 0] (solid line), D q = [0, 0] (dashed line) 
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Figure 24: Bode plot of the closed loop transfer function from a c to a; 
Dq = [-1.99, 0] (solid line), Dq = [0, 0] (dashed line) 


50 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden lor this collection of information is estimated to average 1 hour per response, including the lime lor reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington, VA 22202-4302. and to the Office of Management and Budget, Paperwork Reduction Project (0704-01 S8), Washington, DC 20503. 


1. AGENCY USE ONLY ( Leave blank) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 


February 1995 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


Stable H°° Controller Design for the Longitudinal Dynamics of an Aircraft 


6. AUTHOR(S) 


Hitay Ozbay and Sanjay Garg 


WU-505-62-50 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-9421 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESSES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSOR! NG/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 106847 


11. SUPPLEMENTARY NOTES 


Hitay Ozbay, Ohio State University, Department of Electrical Engineering, 2015 Neil Avenue, Columbus, Ohio 43210; 
Sanjay Gaig, NASA Lewis Research Center. Responsible person, Sanjay Garg, organization code 2550, (216) 433-2355. 


12a. DISTRIBUTION/ AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Categories 63 and 08 

This publication is available from the NASA Center for Aerospace Information, (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 

This report discusses different approaches to stable H°° controller design applied to the problem of augmenting the 
longitudinal dynamics of an aircraft. Stability of the H°° controller is investigated by analyzing the effects of changes in 
the performance index weights, and modifications in the measured outputs. The existence of a stable suboptimal control- 
ler is also investigated. It is shown that this is equivalent to finding a stable controller, whose infinity norm is less than a 
specified bound, for an unstable plant which is determined from parametrization of all H°° controllers. Examples are 
given for a gust alleviation and a command tracking problem. 


14. SUBJECT TERMS 


Multivariable control; Stable controller; Right control 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


15. NUMBER OF PAGES 

53 


16. PRICE CODE 

A04 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 




















